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A scheme for probabilistic entanglement generation between two distant single electron doped 
quantum dots, each placed in a high-Q microcavity, by detecting strong coherent light which has 
interacted dispersively with both subsystems and experienced Faraday rotation due to the spin 
selective trion transitions is discussed. In order to assess the applicability of the scheme for distant 
entanglement generation between atomic qubits proposed by T. D. Ladd et al. [New J. Phys. 8, 184 
(2006)] to two distant quantum dots, one needs to understand the limitations imposed by hyperfine 
interactions of the quantum dot spin with the nuclear spins of the material and by nonidentical 
quantum dots. Feasibility is displayed by calculating the fidelity for Bell state generation analytically 
within an approximate framework. The fidelity is evaluated for a wide range of parameters and 
different pulse lengths, yielding a trade-off between signal and decoherence, as well as a set of 
optimal parameters. Strategies to overcome the effect of non-identical quantum dots on the fidelity 
are examined and the time scales imposed by the nuclear spins are discussed, showing that efficient 
entanglement generation is possible with distant quantum dots. In this context, effects due to light 
hole transitions become important and have to be included. The scheme is discussed for one- as 
well as for two-sided cavities, where one must be careful with reflected light which carries spin 
information. The validity of the approximate method is checked by a more elaborate semiclassical 
simulation which includes trion formation. 
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I. INTRODUCTION 

An electronic spin confined in a semiconductor quan- 
tum dot (QD) is an important candidate for a po tential 
building block of future quantum computergSSEl due to 
the long relaxation and coherence times, which are mea- 
sured to exceed 20 ms and 10 ^s, respectively.^] Entan- 
glement between distant electronic qubits using strong 
coherent light and dispersive interaction has been pro- 
posed to be useful for large distance quantum repeaters.15^ 
In this scheme, direct interactions between the qubits do 
not play a role but the entanglement is achieved by let- 
ting both quantum systems interact with a laser pulse 
which acquires a phase shift conditional on the state of 
the qubit and in turn is measured by homodyne detection 
and the entanglement is distributed over kilometers.!^ 
Thereby, the spin degrees of freedom are projected into 
a maximally entangled state. In Ref. [8] a situation was 
discussed where each electronic qubit is placed in a high- 
Q microcavity for better results. 

In this work, we analyze feasibility of this scheme 
as a laboratory experiment of high-fidelity entangle- 
ment creation using the spin of an excess electron in a 
self-assembled QD in a cavity as qubit and exploiting 
the spin-selective trion transitions which lead to Fara- 
day rotation of the light which has interacted disper- 
sively. Faraday rotation with QDs has been measured 
by Atatiire et alWThe main source of decoherence is due 
to light scattering when it intera cts w ith the QDs, and 
thus we aim to analyse the fidelitjJ^ESl fQj. ggll state gen- 
eration, including the effects of the measurement uncer- 



tainty and of decoherence. As for typical QD parameters 
the saturation of the interaction is rather low, it is very 
helpful to use a simple analytical model by eliminating 
the excited states, which enormously simplifies the anal- 
ysis of the fidelity in terms of all the parameters involved 
and the identification of their optimal values. 

In order to apply the scheme of Ref. [Sjto distant QDs, 
one has to take into account that QDs are less ideal ob- 
jects than atoms and, in general, the two QDs do not have 
equal properties. As entanglement generation relies upon 
indistinguishability of the two cases where the QDs have 
opposite spin, it is important to understand the depen- 
dence on deviating parameters and to work out strategies 
to overcome this limitation. Eventually, we also have to 
consider the valence band structure and, in addition to 
heavy hole transitions, also take into account light hole 
transitions which are further detuned and couple more 
weakly. However, the effects are non-negligible for some 
scenarios. 

As QDs interact with the nuclear spins in the solid 
which let the created entangled states dephase, limiting 
time scales are imposed which will be discussed. Satura- 
tion effects not included in the simple model are checked 
via semiclassical simulations. Primarily, we discuss not 
only one-sided cavities which can be used with acousto- 
optic modulators (AOM) that bring in the laser and then 
allow for the reflected light to go off in a different direc- 
tion but also two-sided cavities which seem simpler in 
the sense that the light can linearly pass them. However, 
the reflected light lost into the environment destroys the 
entanglement and one has to be careful here. 

First, we shortly describe trions and distant entangle- 
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ment generation, motivate the use of a cavity, and de- 
scribe the time scales imposed by the nuclear spins in 
Sec. [n] The systems Hamiltonian, the expansion, and 
approximations made in order to get a simpler model, 
as well as the resulting dynamics, are discussed for sev- 
eral pulse lengths in Sec. |III| followed by the discussion 
and evaluation of the fidelity with equal parameters at 
each QD-cavity system. The problems when using two- 
sided cavities and implications are presented thereafter. 
In Sec. |IV| we go into the issue of nonidentical QDs 
followed by a discussion of light hole transitions and its 
effect on the scheme in Sec. |V] Finally, we test our model 
by semiclassical simulations before we draw conclusions 
in Sec. |VII| The appendices are concerned with details 
about the expression for the fidelity and details related 
to nonidentical QDs. 



II. SCHEME AND ITS LIMITATIONS 

When a gate voltage is applied such that exactly one 
electron tunnels in, a QD has two ground states with 
spin ±i (I t), I D) and thus also two possible excitations 
called trions, consisting of a of hole and two electrons 
with antiparallel spin in the lowest lying conduction band 
state (I tijll)! I TijJJ-)) with spin ±|, when only heavy 
hole transitions are considered. According to the optical 
selection rules, | is excited only with left circular 

(or plus) polarized light, and | tii-IJ-) only with right (or 
minus) circular polarized light (strong trion transitions). 
The net effect of spinflips induced by heavy-light hole 
mixing, leading to nonspin preserving (weak trion) tran- 
sitions and coherently coupling the ground states, is as 
small as that due to nuclear spins at an externally applied 
magnetic field B^xt = 1 T.l^ 

Thus, effectively we can treat the QD as a four-level 
system, where light of a definite circular polarization only 
sees a two-level system. The interaction with a highly de- 
tuned field is mainly dispersive and can be calculated by 
putting the susceptibility x of a nonabsorptive mediunPil 
in the limit of large detuning into the slowly varying en- 
velope approximation (SVEA) equation. The phase shift 
acquired after interaction with the field of a laser pulse 
of length L, with cross section Al and frequency wl, can 
be expressed a^^U 



Hamiltonian for the laser light and QD spins. 
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is the 



where Aco = — wr, is the detuning, T — 5- 

radiative decay rate, fj,eg is the dipole matrix element 
of the transition with frequency 1^, and tJo = is the 
scattering cross section of the QD with A = ^ and 
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When dissipation is neglected, we have an effective 
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at each subsystem x = A^B, where q = 1 stands for 
spin up or plus polarization and q — Q for spin down or 
minus polarization, depending on the context, gq denotes 
the ground state of spin q. When using a;-polarized light 
as input {Mx\Q)y = l7|> + l^>-), we get for J,, = 



J, interaction time t 



and 6* ^ 1 a rotation of the 



polarization plane (Faraday rotation represented by 



UaUb^[\9o)a + \gi)A\ ^ [\go)B - \gi)B] |0), (3) 

-> ^[E \9g)A\gg)B\i-ir0a,^)y + V2\^-my] , 

where we discard the x-polarized component because it 
does not matter. I^-") = -^[\go) A\gi) b - \9i)A\go)B] is 
the Bell singlet state. Experimentally, the initial state of 
the QDs in Eq. (3} can be achieved by spin-flip Raman 
transit ion applied to a spin state prepared with high 
fidelity by spin pumping.S^ 

We will show now by estimated conditions, that over- 
coming photon scattering requires the need of a cavity, 
similar to Ref. [13] where the case of a quantum non- 
demolition measurement using Faraday rotation was dis- 
cussed. We note here that their conditions were certainly 
too strict as scattering is not harmful in a readout experi- 
ment. Decoherence caused by the decay of the trions with 
a linewidth of maximally T — 0.002 meV leads to elastic 
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{Rayleigh) scattering at rate 

is the light-matter coupling) at the emitter,iJl'thereby re- 
vealing the spin state. Therefore, the number of scatter- 
ing processes should be kept small, basically less than 1, 
while the signal-to-noise ratio (SNR) from Eq. ([3} should 
be greater than 1. When a one-sided cavity is used for 
each QD, having the effect of enhancing the coupling by 
the finesse factor = where Trt = is the round- 
trip time, is the lengtn, A^ is the area, and k is the 
decay rate of the cavity, we aim to fulfill the following 
conditions: 



SNR = 



'^scatt 



> 1 



J- CTo Taj, 
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For a cavity with T ~ 10*^, there is a regime due to the, 
respectively, linear and quadratic dependences of SNR 
and Uscatt on , in sharp contrast to the case without a 
cavity. Elimination of aj^ yields the necessary condition. 



9'> 



(5) 



which corresponds to an intermediate coupling regime. 
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Compared to atoms, QDs are certainly less ideal ob- 
jects: the transition frequency and the strength of the 
light-matter interaction of two self-assembled QDs will 
in general never be identical as they cannot be controlled 
in the growing process and the QDs are chosen out of 
many randomly distributed samples. As our parameter 
space is quite large, analyzing strategies to get around 
this problem are helped a lot by using a simpler model 
than solving differential equations for each set of param- 
eters. 

For the preservation of entanglement between distant 
QDs, the limiting time scales are due to hyperfine cou- 
pling with the mesoscopic bath of the nuclear spins of 
the lattice. For an externally applied magnetic field of 
the size of iJe^t ~ 1 T, spin flips are largely suppressed.^! 
The effective magnetic field of the nuclei orthogonal to 
Bext can be eliminated by means of a rotating wave 
approximation.GSIln other words, spin flips are prevented 
by energy conservation. As the nuclear spin correla- 
tion time (~ 1 ms) is large compared to the timescale 
for entanglement generation, we may treat them in the 
quasistatic approximation assuming a constant nuclear 
("Overhauser") field of B„„c = 15 mT for InAs/GaAs- 
QDs,l^ which is different for each experimental run and 
may be treated as a classical Gaussian distributed ran- 
dom variable.li^ Thus, as singlet and triplet zero (sin- 
glet with a relative plus sign) get mixed in an unknown 
way due to the different B:^^^ and B^^^ at each dot, 
the entangled state completely dephases at a timescale of 
^2* = {'leBnuc)~^ 1-26 ns. However, by applying spin 
echo, which should be uncomplicated when using electric- 
dipole-induced spin resonance,!^ the singlet rephases af- 
ter twice the time interval between the preparation of the 
initial state from Eq. ^ and the spin flips: At any given 
time the two QD spins are in an unknown superposition 
of singlet and triplet. When the entangled state is going 
to be used for some task at a specified time, spin echo is 
used to ensure that the state has rephased into a singlet. 
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Hjc 



cav — bath 



HqD-bath 



The cavity mode operators (Xq with energy luq are cou- 
pled with coupling constant ^/k to Markovian reservoirs 



Clearly, one must keep track also of the phases when the 
two QDs differ in parameters and due to Zeeman split- 
ting by Beict- Spin echo signals, as has been measured, 
e.g., by Petta et. aZ.,E3 decay due to the variation of the 
nuclear spins at longer time scales i.e. spin coherence is 
lost irreversibly into the environment at r2 ~ 10 

The time scale on which entanglement can be gener- 
ated is now determined by the time during which the 
initial state can be prepared (tprep < 1 ns)^ by the prop- 
agation time {tprop ~ ns), and by the pulse length which 
we will determine below. 



III. FIDELITY WITHIN AN APPROXIMATE 
MODEL 

Our strategy is to first expand the Hamiltonian 
for a cavity containing a QD spin (four-level system) 
and eliminate the upper levels. From this strongly 
simplified Hamiltonian, we derive a Markovian master 
equation and, since the expansion implies discarding all 
anharmonic terms, treat the light classically. 



A. Hamiltonian 

The Hamiltonian for a single one-sided cavity contain- 
ing a QD with one excess electron and driven by a laser 
pulse is obtained by making the typical approximations 
which are common in quantum optics for a system with 
several inputs and outputs and a microscopic description 
of system and bath,'^^ 

H = Hq+ Hcav-bath + HqD-bath + Hjc , (6) 

with 



(7) 



described by continuum operators 6g(w). We denote the 



q=0,l 9=0,1 9=0,1"^ 

9=04 '' 

Y 9{^t^1 + ^q"-q) ' 
9=0,1 

Y j «'^(«g^9(^) - "9^gM) ' 



9=0,1 



9=0,1 
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four states of the QD as |(7g=o,i) for the ground states 
and 1 6^=0,1 ) for the excited states (trions), separated 
by an energy v. Including the Zeeman shifts we have 
transition frequencies Vq — v — {—lY{gh — ge)lJ-BBext + 
{ — l)'^gefJ-BBnuc, where the electron and hole g factors 
are —0.6 and 1.8, respectively, and the nuclear spins only 
couple to the electrons. As Bnuc ^ B^r^t, its effect on 
the phase shift is negligible. The operators describing 
this four-level system are S+ = \eq){gq\, = \gq){eq\ 
and = \eq){eq\ — \gq){gq\- These are coupled to reser- 
voir operators Cg(a;) with coupling constant -s/T and to 
the cavity fields within the Jaynes-Cummings modePs] 
with g. 



B. Expansion 

A systematic expansion for large detuning is achieved 
by applying a Schrieffer- Wolff transformatioipSlto Eq. ([6| 



detuning between the driving field and the trion transi- 
tions, the population in the excited state is very low and 
we neglect that term. This is the approximation which 
renders the transformed Hamiltonian particularly simple 
because the excited states can be completely eliminated. 
Rayleigh scattering is described by 

(12) 

duj{alcq{uj) - agcliuj)) , 



Huayleigh — [AjHqo^bath 
q=0,l 



2tT AuJa ' 



and provides the main decoherence process, as will be 
discussed in Sec. Jll 

The master equation for the density operator p in an 
interaction picture with respect to Hq is obtained after 
elimination of the excited states by making the Born- 
Markov approximation, common in quantum optics,'22l 
and by discarding fast rotating terms proportional to 
g±iAc^,t^ such as those where Hcav-bath and HRayieigh 
are mixed, 



H+[A,H] + l[A, [A,Ho]] + 0{^] 



Ttq a'qf 



(8) 

^"l and detuning 

Wo (or Auj :— V — LdQ without B field). Typ- 
ically, the coupling constant g is at least 1 order of mag- 
nitude smaller than the detuning and the expansion is 



with A = J2q=0,l^S^t^, 

Aujq := lyq 



an excellent approximation, provided that 



< 1. 



The interaction term between QDs and cavity light fields 
is transformed away in first order due to [A, Hq] = —Hjc 
and the ideal interaction, i.e., Faraday rotation, is con- 
tained in 



(9) 



2 ^ Aw, 

9=0,1 ' 



\Pq + tl{l + 2alaq)] , 



consisting of Lamb and Stark shifts. We defined projec- 
tors onto the subspaces of given spin q Pq = \eq){eq \ + 
\gq){gq\- The resulting Hamiltonian is then 



H = Ho + H 



FR + Hcav-bath + Hgo-bath (10) 
+HpurceU + Huayleigh + ^{ ^^2 ) ' 

The last two contributions of H lead to additional decay 
of the cavity fields and the trions via the interaction. 



H Purcell ■ Hcav-bath 

/ K 5 



(11) 



9=0,1 
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du;{t+bq{Lo)-t-bUco)) 



describes the Purcell effeclPil and leads to driving of the 
trion transitions because of the coherent excitation of the 
reservoir modeled by the bq{uj) fields. Due to the large 



dp 
di 



„ , Awo 



alaq\gq){gq\,p] 



(13) 



9=0,1 



[dq 



9=0,1 

9=0,1 

+ (LqpL\-\{LlLq,p}) , 

9=0,1 

where the first term accounts for the Stark shift, the 
second and third ones describe a driven damped cavity 
and the last one Rayleigh scattering at rate = 

with Lindblad operator Lq = ^Tj^dq\gq) {gq\, i.e., spin- 
dependent light scattering which eventually leads to the 
decay of the coherences. The driving field required for 
Faraday rotation is an x-polarized driving laser pulse cen- 
tered at to and is given the shape 



(14) 



This means, each cir- 



as well as photon amplitude 

cular polarization carries a number of photons Its 
pulse length is Tp and its central frequency is oj^. This 
dynamics clearly implies a classical evolution of the light 
(for the mean values, {\gq){gq\dldq) = {\9q){9q\){aldq) 
certainly holds since \gq){gq\ is constant in time), i.e., all 
terms that would lead to quantum corrections to the light 
are in higher order in the expansion parameter. 



C. Cavity fields and signal 

Making an ansatz for the light of either circular polar- 
ization in terms of coherent states, the density operator 
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for the entangled atom-cavity system is given as 
/5W = J2 {p9^9M9g){9q\<^\ag{t))g{&g{t)\g (15) 



(a) 



(g)\ag'{t))q>(ag.{t)\q^ 



J2 [p9,'gM\9q){9q' 



where q' denotes the polarization opposite to q. aq{t) 
and aq{t), respectively, stand for the amplitudes of the 
coherent states when the QD is in the interacting and 
noninteracting spin. 

The equations of motion for aq{t) and aq{t) are given 
as the derivative of the cavity field mean value (a) = 
Tr{a/5} for the cases Pg^,g^,{^) = 1 and Pg^g^iO) — 1, 
respectively (note that these density matrix elements are 
constant in time as we eliminated the excited states) us- 



ing the master equation [Eq. (13l]. Defining ^^S{t) :— 
cxq{t) for an empty cavity, the cavity field has, when 
driven on resonance, the shape (here and in the follow- 
ing, we will always neglect the damping of the light due 
to Rayleigh scattering as in the regimes of interest this 
effect is negligible) 

« ^e-'"»*^;H(i) , (16) 

which is the solution for a driven harmonic oscillator and 
the approximation holds in steady state (rp ^ i). Going 
on to the case where the cavity field interacts with a QD, 
^^Sq{t) := aq{t) is given on resonance by 

Sq{t) = Vf^e-'^"' rdi'e(-^(»'/^"')+("/2»(*'-*)F™(i') 
Jo 



-iujQt 



(17) 



The output field is related to the cavity field and the 
input field by means of the relation &ouT(i) = y/^aq(t) — 
bj^{t) (Ref. [T8| and is a continuum field operatoil^ with 
the structure 



(18) 



whereby FouT(i) — W^Sq{t) — F-i^{t)\ is the pulse shape 
with ^'^^dt\Fo-uT{i)\^ — 1- The corresponding eigen- 
states are the coherent states. 



a? 



l,,,{t)) = exp / dtlal^^Mb^oUt) - H.c] |0), , 



(19) 

where the integral is over the real axes. For the case of 
dispersive interaction, we have 



V2 



(20) 
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FIG. 1: \S{t)\ (empty cavity) for k = 0.05 meV (gray solid 
line) and output field [Eq. ( 18 l] in units, where k = 1 (light 
gray solid line) compared to the steady state estimates (black 
solid line) for pulse lengths of (a) 1 ns, (b) 100 ps, and (c) 
10 ps. The dashed lines illustrate the phase shift in terms 

2 

of i^^- The parameters are g = 0.15 meV, Aoj = 5 meV, 
and QiN = 8. For rp = 10 ps the rapid change of the phase 
corresponds to a sign change. 



which corresponds to the Faraday rotation to be mea- 
sured and is equal to the estimation from Sec. |ll] In 
order to check the validity of the steady state assump- 
tion, we calculate Eq. ([l6| by solving the differential 
equation fulfilled by S{t) for several pulse lengths with 
results plotted in Fig. [Tiand compare to the steady state 
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curves. The phase shift is plotted in terms of corre- 
sponding to the phase shifts for large detuning. For short 
pulses, the shape gets deformed and translated and we 
encounter nonconstant arguments for short pulse length 
(tp < 1 ns). However, it is always possible to replace 



the nonconstant phase 9q{t) in Eq. (19 1 by a mean phase 
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defined as 



= / dt|F„V(t)l 



2gie,(t) 



(21) 



By carrying out the integral in Eq. (19 1, we can equiva- 
lently denote the output states as |y=e*^'')g, i.e., by the 
photons carried by the pulse. As can be seen in Fig. [l] 
for pulse lengths not shorter than Tp — 100 ps, the steady 
state approximation is quite good. Also, the output pulse 
shape does not depend on q and we discard that index. 
This allows a simple decomposition into linear polarized 
components as 



Using pulses with Tp '-^ 100 ps, one may create en- 
tanglement on a time scale of ~ 10 ns: The interaction 
of the light with one cavity lasts for about ~ 1 ns, the 
light travels typically 1 ns between the cavities, and after 
the interaction, a spin-echo pulse can be applied which 
rephases the desired Bell state after twice the interaction 
time. 



D. Decay of coherences 



The equation describing the coherence 



^out('') ~ -FouT(i)e ^ ^^(aQ^^-t-z( l)'^aQ^^) 



V2 



(5S„,(t)-H*(-l)'^6^„,(i)) 



(22) 



Balanced homodyne detection offers a means to mea- 
sure the quadrature operators of the electromagnetic 
output field by integrating the field of interest with a 
large state (local oscillator) at a beam splitter and sub- 
tracting the photocurrents produced by the outpuLl^il 
Here we consider a situation where a light pulse of 45° 
linear polarization is split into a y-polarized compo- 
nent, which serves as the (classical) local oscillator state 
Q^Lo(i) = Q^iN^iN(i)e~"^^*5 and an a;-polarized compo- 
nent, which interacts with the QD-cavity system whereby 
a y-polarized component may be acquired. The x quadra- 
ture of the latter is to be measured and thus the x- 
polarized component is removed at a polarizing beam 
splitter. The observable representing the homodyne 
detectoi^^ is then given by 



St ^ Qh dt[{bl^^{t)ya^o{t) + H.c.] 



(23) 
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where is a constant related to the measurement appa- 
ratus and the integral is over the entire pulse duration. 
Thus at the homodyne detector, the density matrix in 



Eq. (15 I, which can equivalently be written in terms of 
output fields instead of cavity fields, is projected into 
eigenstates of x^ut — ^(oout + (flouT)^)- The factor 
in front of x^ut, in particular, J dtFjf^{t)FouT{t), does 
not matter as long as it is not too small like for pulse 
lengths Tp < (see Fig. [T||. Then, clearly the relative 
noise is not solely determmed by the variance of x^, 
but other contributions become important which have 
been neglected due to the large interference between 
signal pulse and local oscillator.i^S We will thus restrict 
ourselves to pulse lengths Tp = 1 ns and 100 ps where 
TpK ~ 70 and 7, respectively. Alternatively, the local 
oscillator pulse could be sent through the same but 
empty cavity structure as the signal pulse, then having 
a very similar shape and improved overlap. 



Paiaoit) = rr(^|.9i)(5o|® |«i(t))i(«i(t)|i (24) 
(g)\ao{t))o{ao{t)\op{t)] 



is according to Eq. ( 13 1 given by 



dt 



1 



(ro^ + rf/2) 



a 



){al/2)\S{t)\'p,,,„{t) 
hs{t)\^Pa^9o{t), (25) 



while the diagonal elements of p{t) are constant. We re- 

placed aq{t) by Uqit) since the terms are already C'(^)- 
This equation describes, besides trivial phases due to Zee- 
man splitting (an additional phase should be added be- 
cause magnetic fields also lead to a relative energy bias 
between the two trion transitions but we will neglect 
that in the following), the decay of coherent superpo- 
sition states caused by Rayleigh scattering. Integration 
yields for the modulus 



\Paigo {Tend) I 



|Psi3o(0)|e~ 
|P9iSQ(0)|e" 



(26) 



taken at Tend sufficiently long after the interaction such 
that no scattering occurs anymore. <i> = ^^^dt\S{t)\^ 
is the pulse area, describing how many photons couple 
into the cavity. The decay in steady state, when accord- 
ing to Eq. (16l <i> = ^, is determined by the estimated 

number of scattered photons n^catt — '^a^ which is 
equivalent to the expression Eq. ^ in Sec. [n[ 

For two one-sided cavities, each containing a QD, the 
decay is determined by the sum of all contributions from 
each cavity x with spin q, 



\PgoQ 



, (Tend) I - \Pgo 



,(0)K 



(27) 



In this density matrix element, which is — ^ for a Bell 
singlet and for a product state, the first and third in- 
deces refer to the first subsystem and the others to the 
second one. 
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E. Fidelity 

When the second cavity is driven by the output from 
the first cavity, the signal has y components given by 

Mil = sin (6*^1 +6*51))^ , 

Moo = ^ sin (6*^0 + dBo))y , 



\d 
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-(sin 6*^1 - sin6'Bo))y , 



Moi = -Y~{sineAQ ~ sin6'Bi))y 



(28) 



We refer to the amplitudes du (both spins up), doo (both 
spins down), dio (A: spin up and B: spin down), and doi 
(A: spin down and B: spin up) as distinguishabilitiesW 
They are the centers of the corresponding probability dis- 
tributions 



^qiq-ii.-^) (-''OUT I J/ 1*^9192)1/ 



(29) 



How close to zero dig and dgi will be determined by how 
well the QD cavity parameters for the two subsystems 
match. c?oo and du are on the order of 1 and mark the 
unwanted situation of parallel spins. The phase shift at 
QD X and for polarization q is given approximately by 
(^xa — — 7^ — ■ The more complete expressions which will 
be used in the following are given in Appendix [X| 

After measurement, the density operator for the QDs 
and output fields depends on the outcome of the mea- 
surement (x) and is given as 

/5(a;,i) = {p9,■,gu■^9,^a,2^t) ■\9q19q2) 

91,92,93,94=0,1 

X (593594 |Gqi,2(a;)G,354 (2;)) ■ (30) 
The fidelity for a singlet is defined as 



F = 



f-^^Jx{^-\p{x,T,^a)\^-) 
TrQo Jll dxp{x,Tend) 



(31) 



normalized by the success probability^ The evaluated ex- 
pressions are given in Appendix |B] Obviously, a mea- 
surement window Xc must be chosen to account for the 
overlaps of the Gaussian peaks, which introduce an un- 
certainty in the measurement result, by defining an inter- 
val around x = within which the measurement outcome 
is accepted and outside of which it is discarded. The in- 



tegrals over the Gaussian functions in Eq. (30 1 are given 



by error functions and the diagonal density matrix ele- 
ments of the initial product state Eq. (|3| are constant in 
time at i and Pg.gomgii^) = Pao9i;9m)fi) = ~h 



F. Results for identical quantum dots and cavities 



Equation ( 31 1 is now evaluated and displayed in Fig. 
|2]for the same parameters at each subsystem, chosen to 



he g = 0.15 meV, and k = 0.05 meV, which seems ex- 
perimentally realistic and means a high-fidelity regime 
[see also Fig. |3](b)], while and Aui are varied. The 
range for possible detunings is between 1 and 10 meV, 
as smaller detunings would lead to significant electronic 
excitations, whereas for larger detuning, one drives un- 
wanted remote (light-hole) transitions and we start to 
lose the polarization-spin correspondence, as will be dis- 
cussed in Sec. |V] A measurement window of Xc = 0.3 
leaves a success probability of ^ 25% in the region of 
interest. A magnetic field of = 1 T is used but this 
does not significantly change the results. A detuning of 
Alu < 2 meV is obviously not a good choice as the fidelity 
is bad and becomes also strongly S-field dependent. 
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FIG. 2: (Color online) (a) Fidelity for identical parameters at 
each subsystem g — 0.15 meV and k — 0.05 meV versus am 
and Au>. A trade-off between decoherence and signal gives 
rise to an optimal set of parameters, (b) Success probability 
Ps corresponding to the fidelity in (a), never lower than 25%. 



A trade-off between the signal strength and decoher- 
ence, similar as in Ref. [5] where the situation of photon 
losses for far distant qubits was discussed, compared to 
photon losses due to scattering here, becomes obvious: 
At low signal, decoherence is unimportant but the uncer- 
tainty of the measurement does not allow us to project 
into an entangled state with high probability, whereas at 
higher signal intensity, decoherence via Rayleigh scatter- 
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ing becomes crucial. To establish a link to the estimates 
from Eq. ([i}, we note that the distinguishability dqq cor- 
responds to the SNR and plot it together with the num- 
ber of scattered photons and the resulting fidelity in Fig. 
|3]^a). The criteria SNR > 1 and riscatt < 1 used at the 
beginning were obviously quite good in terms of roughly 
revealing the regime of high fidelity. The necessary con- 
dition for the operation regime estimated in Sec. |ll]to be 
g2 ^ kT pg^jj jjQ^ tested, as shown in Fig. 3[b), where 
we calculated the maximum possible fidelity for varying 

2 

^ (with r fixed at 0.002 meV) which can be compared 
to the estimated condition from Eq. ([5|, yielding at the 
boarder of the inequality F ^ 0.7. For the parameters 

2 

g = 0.15 meV and k = 0.05 meV we have ^ = 225, thus 
they belong to the high-fidelity ^ 0.99 regime. 




0.25 0.5 0.75 1 1.25 1.5 1.75 2 
Xc (arb . units ) 



FIG. 4: Fidelity (solid lines) and success probability (broken 
lines) versus measurement window for various detunings of 
2 meV, 4 meV, 6 meV, and 10 meV corresponding to the 
gray level changing from dark to light gray. 
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FIG. 3: (Color online) (a) Comparison of SNR and total 
number of scattered photons from Sec. [lljand the fidelity for 
Au — 5 meV versus Oin- (b) Maximal possible fidelity for 

2 

different ratios of ^ logarithmically scaled. At the border of 
the inequality Eq. p|, we have F ~ 0.7. 

In more detail, the connection between success prob- 
ability and fidelity is displayed in Fig. |4] showing how 
much fidelity we lose if we require maximal success prob- 
ability. 



G. Discussion of two-sided cavities 

Coupling one-sided cavities requires the use of AOMs 
and thus it seems appealing to use two-sided cavities 
and directly send the transmitted light to another cavity. 



However, we must cope with the fact that the reflected 
light of a double two-sided cavity structure containing 
QDs carries as much information about the spin state as 
the transmitted light. Additionally, there might also be 
internal reflections between the cavities when the laser is 
not perfectly on resonance. A general discussion of this 
case can be found in Ref. |2ll 

We consider a situation with left-incident light. The 
formulas derived for one-sided cavities are adapted by 
considering that the leaking of the light from the cavity 
is twice as high (assuming identical mirrors), while the 
coupling constant ^/k of the output remains the same. 
This decreases the number of photons in the cavity by 
a factor of four. The phase shift of the output in trans- 
mission is the same as that of the cavity field and thus 
4 times smaller as compared to Eq. (20 1. If there are 



reflections, we also have to consider that the signal is de- 
teriorated due to reflected light by the pulse area of the 
second cavity, i.e., by a factor of ^/^bi^b, which comes 
in as the normalization factor of the output [analogous 
to Eq. (18 1 for the case of one-sided cavities]. 

Entanglement generation relies upon detecting the out- 
put of the second cavity where the light carries informa- 
tion about both subsystems. As a two-sided cavity has 
outputs with spin information at any mirror, entangle- 
ment is destroyed when the reflected light can in prin- 
ciple be detected. Due to the boundary conditions at 
one cavity mirror (cf. Sec. IIICl, on resonance (with 



the empty cavity) reflection for light with the polariza- 
tion corresponding to the active spin state and thus the 
sign of the y-polarized component depends on spin orien- 
tation (but the x-polarized component does not). With 

2 

dxq '■— ^f^]^^ — the y-polarized output of the double- 
cavity system at the driven (left) mirror, i.e., the y- 
component of the reflected light, is then approximately, 
for the two spin configurations of interest l^i^o) ^ind 
15051), 



yL 



IdAlFj^it) — dBoPm{t — '2tprop))yL , (32) 

I - dAoFjt^{t) + dsii^iN 
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respectively. Thus, for equal QDs, we observe two subse- 
quent pulses of opposite amplitude, the order of which 
depends on the spin states. Tracing over this degree 
of freedom leads to additional decay of the fidelity due 



to the decay of {ps^'io^^^^pg 



lyL 



with details given in 



Appendix |B] If the pulse length is short enough such 
that the reflected pulses do not overlap, the effect has 
its maximum and the fidelity practically decays com- 
pletely. For longer pulse lengths, the overlap integral 
Ioi{tp) ■= J Fj^{t)Fjf^{t — 2tprop)dt approaches unity [see 
Fig. [sja)]. The dependence of the fidelity on the pulse 
length is shown in Fig. for small (2 meV) and large 
(10 meV) detunings with a success probability > 40%.^ 



(a) 



I 

0.8 
I 0.6 
iO.4 




FIG. 5: (a) Overlap integral loi versus pulse length, (b) Maxi- 
mal possible fidelity for a two-sided cavity versus pulse length 
Tp for Au) — 2 meV (gray) and Auj = 10 meV (black) for 
Xc ~ 0.7, yielding Psucc > 40%. Smaller Xc do not apprecia- 
bly increase the fidelity. 



IV. NONIDENTICAL QUANTUM DOTS 

In experiments, the ideal results from the last section 
will not apply since it is not very likely to find two self- 
assembled QDs of same frequency and same g. Here, we 
show how to overcome this problem by several different 
strategies depending on the difference in trion transition 
energies va and vb of the two quantum dots. 

For large detuning of the QDs | va — i'b \, one can tune 
the laser symmetrically in between the QD resonances 
instead of redshifting in order to balance the detuning 
for the two subsystems. Then, one produces the Bell 
triplet state |(/)+) — ^(|00) + |11)) instead of the singlet 
state since the sign of the phase shift depends on the 
sign of the detuning. Thus, the phase shift of the first 



cavity differs by a sign from that of the second cavity and 
the role of |01), |10) and |00), |11) is reversed (see also 
Appendix [A]|. With the help of the formalism introduced 
in the preceding sections with the simple analytic formula 
for the fidelity, it is now easy to analyse the dependence 
of the fidelity on varying Auj a and Aujb by finding the 
optimal parameters numerically, as shown in Fig. |6] The 
resulting fidelity becomes bad when |z^^ — i'bI < 2 meV 
due to Rayleigh scattering, as shown in Fig. |6](blue line). 

In order to have good fidelity when the QDs differ 
not too much, it is advantageous to redshift the cav- 
ity frequency such that it is detuned from the QD with 
the lower transition frequency by the maximal allowed 
value, which we estimated in Sec. |IIIF| to be about 
Auj = 10 meV. In this case, the relative difference of the 
phase shifts at each QD is minimized, while light-hole 
mixing effects remain small. By optimizing the photon 
number, we get high fidelity for {i^a — i^b\ < 2 meV, as 
shown by the red line in Fig. |6] 

For small detunings it may also pay to slightly detune 
the cavity with the "better" (higher transition energy) 
QD from the laser such that the phase shift at this QD 
becomes of the same size as that one at the other, "worse" 
QD (see Appendix |A] for details). We find now optimal 
and laser-cavity detuning at one of the subsystems 
(here A) 5uj — ujl — wqa by optimizing the fidelity (green 
line). However, by doing this, the signal which distin- 
guishes the entangled state from product states \dio)y 
and \dQi)y (cf. Sec. Illl, which is the sum of both phase 



shifts at either cavity, decreases which lowers the SNR. 
Thus, this method does not work for arbitrary Suj. The 
same effect could be achieved by decreasing the measure- 
ment window Xc from 0.3 to 0.2, which, however, de- 
creases the success probability (violet line). 

Realistically, we have also different QD-cavity coupling 
constants gA and gB and we can also compensate for 
this by either tuning in between for big {i^a — i^bI or 
redshifting otherwise if gA < gB and va < vb- Thereby, 
we compensate for, e.g., a smaller gA with a smaller Auja 
(black lines). If gA > gs and va < vb and thus the 
QDs are even more nonidentical, tuning in between works 
well, but for red shifting a smaller measurement window 
Xc = 1 must be chosen. 

In order to have higher success probability, we choose 
a larger measurement window of a:c = 1 and examine 
the same strategies as before in Fig. [7] While the suc- 
cess probability can now be increased above Pgucc = 0.47, 
simple redshifting (red line) is now less suited to compen- 
sate for different QDs. Instead, it pays now to detune one 
cavity slightly from the laser (green line). Tuning in be- 
tween for different QD transition frequencies or different 
light-matter coupling constants still works well. 



V. LIGHT-HOLE TRANSITIONS 

So far, we considered only the dominant heavy hole 
(trion) transitions which is certainly justified for small 
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FIG. 6: (Color online) (a) Outline of the fidelity versus \va — i^b\ for different strategies. Tuning in between the two QD 
resonances does not work for small \ua — i^b\ (blue line), where redshifting is more appropriate (red line) with maximal detuning 
being 10 meV. One could also detune the cavity containing the QD with larger frequency (smaller detuning when redshifting) 
slightly from the laser (green line) or chose a smaller measurement window of Xc = 0.2 (violet line). Then, one obtains a smaller 
success probability as compared to Xc = 0.3, which was the choice for all other lines. The black lines are for qa ~ 0.14 meV 
with qb = 0.15 meV (solid black line) and qb = 0.14 meV with qa ~ 0.15 meV (broken black line), respectively. The second 
case deviates a lot from the first for red shifting and we should accept a lower success probability {xc — 0.1) yielding better 
fidelity (broken violet line), whereas for tuning in between, the lines are indistinguishable. In (b), the optimal photon number 
is shown which is of course bigger when the laser is detuned from the cavity (green line). When qa Qb, the laser is not 
tuned in symmetrically between the QD transition energies but asymmetrically with Aujasy = , as shown in 

(c) (black lines), together with the optimized cavity-laser detuning 5uj = \ul>o — <^l| at one cavity (green line). All the lines are 
interpolations of the evaluated points. In (d), we plot the success probability. 



detuning. Going to higher detuning or when considering 
nonidentical QDs and tuning in between, as in Sec. |IV| 
requires taking into account the valence band structure, 
i.e., also the light-hole transitions. Effects from other, 
remote transitions are discussed in Ref. [25] and were 
found to be small for a detuning of less than 10 meV. 
The light-hole-associated energy levels are separated by 
at least Aujhl — 10 meV from the top heavy-hole level. 
Moreover, their dipole matrix elements are reduced by 
a factor of ^ with respect to to heavy holes, see, e.g., 

Ref. g6] leading to 5 ^ ^ and T ^ §. The allowed 
transitions lead to a coupling of each spin to both circular 
polarizations, as depicted in Fig. |8] 

The Faraday rotation is thus, for a fixed spin state, 
reduced as now both circular polarizations acquire a fi- 
nite but different phase, such that the relative phase de- 
creases. The Stark shift becomes reduced as 



Applying a Schrieffer- Wolff transformation to the Hamil- 
tonian including the light holes, one gets analogously to 
Sec. |III| Rayleigh scattering contributions for all transi- 
tions involved, but we now have two scattering contribu- 
tions for each polarization. This corresponds to choosing 
the same reservoirs for the transitions which are coupled 
by the same polarization. The Lindblad operatoi^SS for a 
certain circular polarization then reads 



L„ 



l59>(5g 
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-\9q'){9q'\)aq 



(34) 



where Aa> — Aojhl ± Auj for either redshifting or 
blueshifting, respectively. This implies that scattering 
of circular polarized light occurs for both spin states but 
at different rates and thus the scattered photons carry 
less information about the spin state as compared to the 
case without light holes. The actual rate at which coher- 
ences decay [cf. Eq. (27l] is decreased and given by the 
replacement 
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FIG. 7: (Color online) (a) Outline of the fidelity versus A = — i^bI for different strategies with same encoding as in Fig. [6] 
A larger measurement window of Xc = 1 was chosen, such that the success probability in (d) increases significantly, with values 
above Paucc = 0.47. In (b) and (c), the calculated optimal photon numbers and detunings, respectively, are shown. 




FIG. 8: (Color online) Level scheme including the light holes 
which have different spins S than the heavy holes. Each spin 
ground state is now addressed by light of both circular polar- 
izations but with different coupling strengths. 



For redshifting, the corrections due to light hole tran- 
sitions are rather small for detunings Auj < 10 meV 
and there is still always a region of high fidelity for 
Auj > 10 meV [see Fig. |9]^a)]. The red region signifies 
where the fidelity is higher than 0.99, while Pgucc > 0.47 
for both the cases of light holes split by Aujhl = 10 meV 
and no light holes. The light yellow region signifies where 
this is true only for the latter case, the darker yellow re- 
gion only for the case with light holes included. Thus, for 
redshifting, we do not have to worry about light holes also 
for bigger detunings. They become more crucial when 
considering nonidentical QDs and tuning in between the 



transition frequencies as in Sec. IV because then one 
QD is blueshifted with respect to the laser and is thus 
close to the light holes. The regime of high fidelity exists 
only for small detuning and photon number, as can be 
seen in Fig. for Psucc > 0.35. Although we could 

find regimes of high fidelity also when light holes are in- 
cluded, it remains the problem that we do not exactly 
know how much they are energetically split from heavy 
holes. Thus, we search for the overlap of the results for 
a pessimistic estimate, Aujhl = 10 meV, and a rather 
optimistic one, Aluhl = 20 meV. In Fig. [lOjwe plot the 
region of high fidelity > 0.99 with a high success prob- 
ability of Psucc > 0.49. Thus, we conclude that a fairly 
high fidelity together with a good success probability can 
be obtained also for the case of different quantum dots 
and the presence of light holes. Experimentally, it should 
be feasible to find two QDs which differ between 2 meV 
and 6 meV and tune in between. The discussions of this 
section are very convincing that QDs represent excellent 
systems for distant entanglement generation. 



VI. SEMICLASSICAL SIMULATION 

We now test the validity of the expansion in Sec. [ll|by 
a semiclassical simulation retaining the excited states, 
similar to Ref. |8l We will consider here for simplicity 
only the situation of identical QDs, Bcxt = 0, and no 
light holes. For the density operator of a single QD, we 
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FIG. 9: (Color online) The colored regions signify where 
F > 0.99. The red region indicates the overlap of the case 
with light holes included {Aujhl ~ 10 meV) and no light 
holes. The light yellow region [left yellow region in (a) and 
in (b)] is only for the latter case and the darker yellow region 
for the first case. In (a), both QDs are redshifted equally and 
Psucc > 0.47 {xc ~ 1), whereas in (b), the laser is tuned in 
between the two QDs and thus one of them is blueshifted, 
restricting the region of high fidelity and Psucc > 0.35 {xc = 
0.6). 
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FIG. 10: The black region signifies where F > 0.99 and the 
gray region F > 0.98 for the overlap region of Aluhl = 
10 meV and Alohl = 20 meV and tuning in between the 
two QDs. The success probability is > 0.49 {xc = 1.3). 



make a similar semiclassical ansatz as in Eq. (15 1, but 
include also the excited states. 



(e<jl + PcgC,{t)\eq){eq\) ® \&g{t))g{&q{t)\q 



9=0,1 



+Ps„/e,(i)|eq)(ffq'| + Pe„, S, (^) ISq ) (e,' | + Pe , e, (t ) | C,) 



® \&qit))q{aq{i)\q ® ! W ) W W I 



(36) 



This ansatz is based on the lowest order approximation 
neglecting any quantum correlations which, according to 
our findings in Sec. |III| and Ref. [8] is a good approxi- 
mation in the low saturation regime. First, we transform 
the Hamiltonian from Eq. ^ to an interaction picture 
with respect to 



\^uJaalag+'^\eg){eg\ + {'^-uJa)\gq){gq\ 

9=0,1 



(37) 



which amounts to the replacements 



a„ 



^ ^q 1 

yielding for the trivial part of the Hamiltonian 



= - ^ ^u,\gq){gq\ 

9=0,1 



(38) 



(39) 



The equations of motion that determine aq{t), i.e., the 
output fields for a definite spin state g, are 



aq{t) 



«.9 Pe,5,W(ag(t))* - (pe,g,(t))*<5,(t) 



(40) 



r ■ 



a 



it) 



--ag{t) + y/^-^F,^{t) - igpe,g,{t) 
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we plot a typical so- 
100 ps, a = 4, and 



where Pg,g,(0) = 1. In Fig. 
lution to this equation for r, 
Aw — 2 meV. The actual value for the phase shift of the 
light is given as an average as discussed in Sec. |III C| and 
the fast increase at short times displays how the steady 
state of the cavity is reached. The phase shifts are com- 
pared in Fig, 
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^a) to the approximate result for different 
pulse lengths, saturating for large values of our expansion 
parameter. However, in the regime Aw = 2 — 10 meV the 
deviation is not more than 10%. 



A test of Eq. (27 1 requires, for a fixed g, solving two 



additional equations for the simplest case of driving with 
circular polarized light (here, g = 0). In this case, we 
have a two- level system {\go), |eo)) driven by nonreso- 
nant light. We are interested in the coherence with the 
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FIG. 11: (Color online) The red line displays a solution of Eq. 



( [iOj l for the phase shift ^^(t) defined in Eq. d20|. The solid 
black line corresponds to Pegeg(i) from Eq. (41 1 showing that 
the population of the excited state is low, whereas the dashed 
line corresponds to pg^e^it)- These amplitudes adiabatically 
follow the cavity field (gray line). Parameters as in text. 



other (ground state) level \gi) which does not couple for 
this polarization. The coherence between the two ground 
states, i.e.. 



,(t), is coupled to pg^eo{t) via 
= igaa{t)pg^eo{t) + i^^Pgtgo{t) ^ (41) 

= w[ao{t)]* Pg^goit) - T^Pgieoit) > (42) 



PgigoW 

Pgieoit) 

where ao{t) :— ^[aQ{t) + aQ{t)] is the solution to Eq. 
(40 1 for Pg(,go(0) — Pgigi(O) — \ accordlug to the situa- 
tion of an initial equal superposition of both spin ground 
states. The approximate and the semiclassical results for 
the damping are shown in Fig. [l2jb) for various detun- 
ings and pulse lengths. 

For the double cavity system, the phase shifts from 
each cavity are added and the total decay of the entan- 
glement coherence is determined by the contributions 
from each transition. Although there is, particularly for 
low Aw, an overestimation of scattering by Eq. (27 1, 



which is due the fact that a smaller amount of light 
couples into the cavity due to an intensity dependent 
Stark shift, we encounter that for the final fidelity, the 
two methods are yielding practically identical results 
[see Fig. [T2jc)]. In Fig. 12 d), the case of tuning the 
laser in between the two QD resonances is shown (see 
also Appendix [X] for the differences between the two 
cases). The saturation of the signal, lowering the fidelity, 
is compensated by a lower decay of the coherence. These 
results suggest that using the approximate model, which 
considerably simplifies practical calculations, is well 
justified. 

Nonclassical effects have not been indnHpH pn far. Go- 
there is a 



IIIB 



ing to C(^^) in the expansion of Sec. 
term describing nonlinear dephasing (see, e.g.', Ref. 



which is 2 orders of magnitude smaller than the linear 
Stark shift. We simulated quadrature squeezing due to 
nonlinear dephasing and found practically no effect for 
parameters where the fidelity is high in Fig. 12 (~ 0.1%). 
This corresponds to the findings in Ref. l^^For short 
pulses, low detuning and high photon number squeezing 
may occur. In more detail, for Aw = 1 meV and a^^ > 4, 
there is squeezing on the order of a percent, whereas for 
higher detuning, it is negligible. Principally, squeezing 
would not harm entanglement creation unless one could 
learn about the spin state at one QD from the amount 
of squeezing observed. However, for our conclusions non- 
classical behavior of the light is not relevant. 



VII. SUMMARY AND CONCLUSIONS 

We have shown that entanglement with high fidelity 
(> 0.99) and success probability (> 0.49) is possible be- 
tween distant quantum dot spins using cavity QED and 
coherent light bus modes. Nuclear spins should not mat- 
ter when spin-echo techniques are used, except for the 
T2 time they impose. A simple analytic model based on 
the elimination of the excited states was used, largely 
simplifying the determination of the optimal parameters 
which are necessary to achieve high fidelity. QDs are gen- 
erally nonidentical in terms of transition frequency and 
light-matter coupling constant, but we have shown that 
there exist strategies which allow to largely compensate 
for this. Taking into account QD-specific effects, such as 
light-hole transitions which mainly become important if 
one tunes the laser in between the two QD frequencies, 
we demonstrate that there exist regimes where all the 
requirements are fulfilled. Thus, effects from the non- 
identical nature of QDs can be overcome. We mainly 
focused not only on one-sided cavities as they are prob- 
ably most suited but also discussed two-sided cavities 
where long enough pulse lengths have to be used, such 
that the spin states are not revealed by the scattered 
light. We tested the simple model against a semiclassical 
simulation which accounts for excited state population 
for several pulse lengths and found it to be valid in the 
regime of interest. 
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FIG. 12: (Color online) Results of a semiclassical simulation, (a) Phase shift from a simulation in terms of the approximate 
phase shift arg[l/{^ ^*A!ir)J "^^I'sus Auj for pulse lengths 1 ns (stars) and 100 ps (circles), and for Qin = 3 (highest), 6, 9, 12, 
and 15 (lowest). For low detuning, high photon number, and short pulse length, we see a strong saturation of the phase 
shift, (b) Damping of coherence Pgogi;gigo{Tend) at one transition with same encoding as in (a) for Au = 1 meV (lowest), 
2 meV, 4 meV, 6 meV, 8 meV, and 10 meV (highest). For low detuning (and thus large Stark shift) the scattering is 
lower compared to the approximate results from Eq. (27l (light blue line) due to reflection induced by the Stark shift, (c) 
Fidelity {xc = 0.3) for different pulse lengths (symbols) compared to the model from Sec. Ill (lines) for Alu — 1 meV (lowest). 



2 meV, 4 meV, 6 meV, 8 meV, and 10 meV (from left to right). In the high-fldelity regime, the results coincide and the large 
deviations for Au = 1 meV are typical for very big Stark shifts ~ ^ [see Eq. (28 1]. In (d), the case of tuning in between the 
QD resonances is displayed. From left to right the lines correspond to Ii^a — i^/2 — 1 meV, 2 meV, 4 meV, 6 meV, and 
8 meV. The colored lines are guides to the eye. 



APPENDIX A: MORE GENERAL EXPRESSIONS FOR THE DISTINGUISHABILITIES 



We explicitly quote here the output states for the general case including slight detuning of the first cavity from the 
laser field by Suj — ojl — ujqa, as used for nonidentical QDs in Sec. IV The four different y-polarized output states 



analogous to Eq. (28 1 are then \d, 



of the light with circular polarization q, 



2 



Ini(7^ ~ 7o n ))v with the corresponding complex amplitudes 7^ 



92 



Iqq 



2 




- 1 




(Al) 



- 1 



The first two expressions correspond to equal spins, while the others to opposite spins. The notion of abrupt changing 
cavity-laser detunings 6uj found by numerical optimization of the fidelity (see Figs. |6]and|7| can be understood from 
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a 




Scj (meV ) 



0.04 



FIG. 13: (Color online) arg(7io) for B = and Aloai = 8 meV (red), Aloai = 6 meV (green), and Aloai = 4 meV (blue) 
compared to arg(7io) for Alobi ~ 10 meV (black). There are two possible solutions for Su) where the lines intersect such that 
the Faraday rotation for both configurations with opposite spins is equal. 



the requirement that 7^ , = 7' , . As can be seen from Fig. 



13 



for a fixed {i^a — i'bI, there are two points where this is 



true: one for positive Slo and one for negative one. The curves are clearly not symmetric in Slo because of the Stark 
shift and thus it is either more favorable to have negative or positive Sus. As soon as the positive solution becomes 
more favorable there is a sudden change. 

For tuning the laser in between the two QD resonances, the amplitudes are obtained by putting a negative Stark 
shift at the second cavity. Thus, for the configuration with same spins, we have for identical quantum dots zero phase 
shifts and finite one for the other configurations. 



APPENDIX B: EXPRESSIONS FOR THE FIDELITY 



The expression of the fidelity defined in Eq. (31 1 gives evaluated 



F = 



(1/2) /^;^ dx{[Gl0{x)]'^Pg„g^,g,g,{Tend) + [^0 1 ( ] ^ ;g„ g„ (^e»d ) " 26*10 (x) Gq j (x) RB ; g„ (Tg^rf) ] } 

„p \ . { ^{erf [V2(a;, - d^)] + erf[V2{x, + dio)] 
erf[V2{xc - doi)] + erf [^2(2;^ + doi)]} 



g-(l/2)(rfio-<ini)' 

with the success probability 

Psucc (-^c) 



dio + doi + 2xc 

7i 



erf 



dio + doi - 2xc 

7i 



-("fN/2) 



dx{[Gll{x)]'^ Pg^g^.g^g^{Tend) + [GqO ( ^ ) ] ^ ^g^ g^ ; g^, g^ ( Tgn ) 

Pgogo;gigi 

{Tend) + [Gf,l{x)Y 

Pgigwgoga {Tend)} 
\{eri[^/2{xe - dii)] + eri[^/2{xc + ^ii)] + erf [^2(2;^ - doo)] 

o 

+CT:i[V2{xe + doo)] + erf [\/2(a;c - d^)] + erf [\/2(xc + dio)] 

+CT:i[V2{Xe - doi)] + GYi[V2{Xe + doi)]} . 



(B2) 



For the two-sided cavity scenario, the ansatz for the density operator [cf. Eq. ( 15 I for the one-sided case] contains 
the photons transmitted and reflected from the double-cavity system. Using Eq. (32 I for the reflected light when the 



spins are opposite, we obtain after tracing out the reflected light the fidelity by replacing in Eq. (Bl I 

giigigo {Tend)] Rc[pg„ {Tend){l3^°^^\yL\P^^^'')yL] , 

with 



(B3) 



l^pSagi ly^l/JfiSf )g£_ = e~(l/2)Exg d^^-2Iai{-rp){dAidBO+dAodBl)\-dAodAl-dBodBl+Ial(rr)(dAodBO+dBldAl) (B4) 
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here dxq and Ioi{tp) are defined as in Sec. 



IIIG 
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